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Heat transport in rotating convection without Ekman layers 
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Numerical simulation of rotating convection in plane layers with free slip boundaries show that the 
convective flows can be classified according to a quantity constructed from the Reynolds, Prandtl and 
Ekman numbers. Three different flow regimes appear: Laminar flow close to the onset of convection, 
turbulent flow in which the heat flow approaches the heat flow of non-rotating convection, and an 
intermediate regime in which the heat flow scales according to a power law independent of thermal 
diffusivity and kinematic viscosity. 

PACS numbers: 47.27.te, 44.25. +f, 47.32.-y, 91.25.Za 



It is a central problem for many areas of geo- and as- 
trophysics to determine the heat flux through a rotating 
and convecting fluid layer. For example, the heat flux 
through the atmosphere governs weather and climate, 
the heat flux through stellar atmospheres determines stel- 
lar evolution, and the heat flux through planetary cores 
is essential for the generation of the magnetic fields of 
these bodies. A correspondingly large effort has already 
been spent on the problem. Buoyancy is driving the flow 
and can be balanced by either viscous or Coriolis forces, 
or the nonlinear terms in the equations of motion, or 
any combination of these. If the Coriolis force dominates 
the dynamics, a special type of boundary layer appears 
near solid boundaries, the Ekman layers, in which the 
viscous force is balanced by the Coriolis term. In ad- 
dition, the flow in the bulk is organized into columnar 
vortices with their axes aligned with the rotation axis. If 
on the contrary nonlinear advection supersedes the Cori- 
olis term, these columns are broken up and the style of 
flow known from non-rotating convection is approached 
[E S Hi ■ There is an ongoing debate concerning the pa- 
rameters at which the transition between these two flow 
regimes occurs 0, 0, @ > an d we are still lacking reliable 
relations between the heat flux and the control parame- 
ters of the flow that would allow us to extrapolate data 
from laboratory experiments and numerical simulation to 
astrophysical objects. 

Some recent work on rotating convection has focused 
on the Ekman layers. For instance, ref. Q relates the 
Ekman layers to the transition mentioned above. Despite 
the inhibiting effect of rotation on turbulence, the heat 
flux in a rotating flow can exceed that of a non-rotating 
flow at equal Rayleigh number p|. In ref. 0, this phe- 
nomenon is attributed to so called Ekman pumps, a term 
reserved for a certain flow pattern associated with Ekman 
boundary layers 0] . Here we investigate convection with 
free slip boundary conditions. This eliminates Ekman 
layers and one can discern which effect really depends on 
their presence. Free slip boundaries are realized to a good 
approximation in Nature, for example at the surface of 
the oceans or at the top of atmospheric layers. 

Consider a plane layer of thickness d in the z— direction 
and of infinite extent in the x, y— plane. Let the layer be 



filled with fluid of kinematic viscosity v, thermal diffusiv- 
ity k, and thermal expansion coefficient a. Gravitational 
acceleration g is pointing in the negative z— direction and 
the layer is rotating with angular velocity f2 about the 
z— axis. The temperatures of the top and bottom bound- 
aries are fixed at To and To + AT, respectively. These 
two boundaries are assumed to be free slip, whereas pe- 
riodic boundary conditions are applied in the x— and 
y— directions. The equations of evolution are made non- 
dimensional by using d 2 / k, d and AT for units of time, 
length, and temperature, respectively. These equations 
then become within the Boussinesq approximation for the 
dimcnsionless velocity v(r,t) and temperature T(r,t): 



Pr 

d t v + (v-\7)v + 2 z xv = -Vp + Pr V 2 v + Ra Pr Tz 

Ek 

(1) 



V • v = 



d t T + v ■ VT = V 2 T 



(2) 



(3) 



z is the unit vector in z-direction and p collects the pres- 
sure and the centrifugal acceleration. The boundary con- 
ditions require that T(z — 0) = 1, T(z = 1) = 0, and 
that v z = d z v x — d z v y — at both z — and 2=1. 
Three independent dimcnsionless control parameters ap- 
pear: The Rayleigh number Ra, the Ekman number Ek, 
and the Prandtl number Pr. They are defined by: 
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The Reynolds number Re and the Nussclt number Nu 
are an output of the simulations: 



, Nu = -jJ <d z T> dA 
(5) 

The angular brackets denote average over time and the 
integrals extend over the computational volume V for Re 
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FIG. 1: Nu as a function of Ra for Pr = 7 (red symbols and 
continuous line) and 0.7 (blue symbols and dot dashed line), 
and Ek = 2x 10" 2 (diamonds), 2 x 10" 3 (squares), 2 x 1CT 4 
(triangles) and 2 x 10~ 5 (stars). Zero rotation is indicated by 
circles and the power law fits have an exponent of 0.287. 
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FIG. 2: (Nu - l)Ek 1/3 as a function of RePrEk 1/2 for the 
same data and with the same symbols as in fig. [1] The 
dashed lines are power laws with exponents 2 and 2/3. The 
two vertical lines indicate the interval outside of which the fit 
to one of the two power laws is considered satisfactory. 



and over the surface A of either the top or the bottom 
boundary for Nu. 

The equations of motion were solved with the same 
spectral method as used in except that free slip 
boundaries were implemented and that the Coriolis term 
was added and treated implicitly together with the dif- 
fusion terms. Resolutions reached up to 129 Chebychev 
polynomials for the discretization of the z— coordinate 
and 256 x 256 Fourier modes in the (x, y)— plane. The 
periodicity lengths along the x— and y— directions were 
always chosen to be identical. The aspect ratio, defined 
as the ratio of the periodicity length in the (x, y)— plane 
and the layer height, was fixed at 10 for simulations with- 
out rotation. In rotating convection, the typical size of 
flow structures varies considerably as a function of the 
control parameters, so that it is not useful to use a single 
aspect ratio. Instead, the aspect ratio was adjusted for 
each Ek to fit at least 8 columnar vortices along both 
the x— and y— directions at the onset of convection, and 
kept constant as Pr and Ra were varied. 

Fig. [T]shows Nu as a function of Ra for various Ek and 
two different Pr. The case of zero rotation is included 
for comparison. The basic features visible in this figure 
are known from previous experiments and simulations 1. 
SB EH- The onset of convection is delayed by rotation. 
After onset, Nu rises more steeply as a function of Ra 
than in the non-rotating case. Nu does not follow any 
simple power in this range of Ra. For large enough Ra, 
the Nu dependence asymptotes towards the dependence 
valid for zero rotation, which is well approximated by a 
power law in the investigated range of Ra. 

All the different curves in fig. Q] collapse to a single 
curve in most of the parameter range when (Nu-l)Ek 1 / 3 
is plotted as a function of Re Pr Ek 1 / 2 as shown in fig. 



[2] For large values of Re Pr Ek 1 / 2 one finds (Nu — 
VfEk 1 ' 3 oc (RePrEk 1 ' 2 ) 2 / 3 or {Nu- 1) oc (RePr) 2 ' 3 . 
This law is independent of Ek as it should be: At any 
fixed Ek and Pr, the limit of large Re corresponds to 
the situation in which the nonlinear term dominates the 
Coriolis term, so that one has to recover the behavior of 
non-rotating convection, which is of course independent 
of Ek. The data for zero rotation cannot be included in 
fig. because Ek has no finite value in this case, but 
(Nu — 1) oc (Re Pr) 2 / 3 is also found for strictly zero 
rotation. 

Low values of Re Pr Ek 1 / 2 on the other hand cor- 
respond to laminar flows near the onset of convection. 
Forming the dot product of eq. |T]) and v, integrating 
over the whole volume and averaging over time, one finds 

e = (Nu — l)Ra (6) 

where e = y J < (diVj)(diVj) > dV is the adimensional 
average dissipation rate of kinetic energy. In a laminar 
flow, one expects e oc (RePr) 2 /X 2 , where A is a charac- 
teristic length scale of the flow. For Pr > 0.676, con- 
vection starts at a critical Rayleigh number Ra c obeying 
Ra c oc Pfc _4//3 and forms stationary cells of size A c with 
A c oc Pfc 1 / 3 [llj]. Eq. © becomes Ek- 2 ^ 3 Re 2 Pr 2 oc 
(Nu — l)Ra. Close to onset, Ra « Ra c and therefore 
(Nu — 1) oc Re 2 Pr 2 Ek 2 / 3 . This corresponds to the 
left asymptote in fig. [21 Both the left asymptote and 
(Nu — 1) oc (Re Pr ) 2 / 3 become straight lines in a loga- 
rithmic plot of (Nu - VjEk 1 / 3 vs. RePrEk 1 / 2 , which 
explains the simple appearance of fig. [21 

Fig. [5] in summary identifies three regimes of rotat- 
ing convection. Rotating laminar flow characterizes one 
of them, and heat transport behaves the same as in 
non-rotating convection in another. The transition oc- 
curs where the two asymptotes in fig. [5] cross, i.e. at 
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FIG. 3: \ m Ek~ 1/z as a function of RePrEk 1 ' 2 . The symbols 
have the same meaning as in fig. [T] The dashed lines indicate 
power laws with exponents and -1/2. 



Re Pr Ek 1 / 2 — 2. There is a transition interval around 
this point of about one decade in width in which Nu is 
close to neither asymptote. This third regime will receive 
detailed attention below. 

Even though Nu behaves as if there was no rotation for 
Re Pr Ek 1 / 2 > 10 in fig. visualizations of the flow still 
reveal differences. In the rotating case, the flow forms 
columnar vortices extending from one boundary to the 
other, whereas for zero rotation, plumes advected by a 
large scale circulation are observed. Enough visualiza- 
tions of vortices in rotating convection have already ap- 
peared [3, 1 1 so that there is no need to reproduce 
any here. The size of the vortices can be quantified by 
the method already used in |9|: Compute the time aver- 
aged advective heat transport through the plane z = 0.5, 
< v z & >, with 6 = T — 1/2, compute the Fourier trans- 
form of < v z Q >, and plot the spectrum of < v z Q > as a 
function of wavelength A (see Q for detailed formulas). 
The median wavelength A m is extracted from the spec- 
tra, such that the heat advected at wavelengths smaller 
than A m equals the heat advected at larger wavelengths. 
The value of A m /2 matches the diameter of the columnar 
vortices identified visually in the flow field. Fig. [3] shows 
XmEk -1 / 3, as a function of RePrEk 1 / 2 . It is seen that 
A m stays at the onset wavelength A c well into the tran- 
sition interval and decreases at high Re Pr Ek 1 / 2 . This 
decrease follows A m oc (Re Pr)~ x l 2 at fixed Ek, which is 
compatible with experimental data in (T3 |. 

Near the onset of convection, the heat transport is de- 
termined by a balance between buoyancy, Coriolis and 
diffusive terms. For high Re Pr Ek 1 / 2 , the Coriolis term 
is overwhelmed by the nonlinear term in eq. ([T]) so that 
Nu is the same as in turbulent, non-rotating convec- 
tion. Diffusive processes play a role because all heat 
has to cross the thermal boundary layers diffusively. Let 
us assume as a working hypothesis that the heat flow 
in the intermediate regime is governed by a competi- 
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FIG. 4: Nu* as a function of Raf*. The symbols have the 
same meaning as in fig. [T] This figure contains only those 
data points which lie in the interval marked by vertical lines 
in fig. H 



tion between the nonlinear and Coriolis terms, and that 
the constraints imposed by rotation on the flow struc- 
ture control the heat flux, not diffusion in the bound- 
ary layers. The dimensional heat flow Q must then be 
given by an expression independent of k and v. In or- 
der to check this hypothesis, it is convenient to use a 
control parameter independent of k and v. The only 
combination of Ra, Ek and Pr meeting this require- 
ment is Ra* — RaEk 2 /Pr = gaAT/(fl 2 d). An appro- 
priate measure of heat flux independent of k and v is 
Nu* — Nu Ek/ Pr = Q / (pc p ATild), in which p stands 
for the density and c p for the heat capacity. 

It is useful to replace Ra* by the flux Rayleigh num- 
ber Raf* given by Raf* — Ra*Nu* — (gaQ) / (pc p fl 3 d 2 ) . 
This combination is strictly speaking a control parame- 
ter only when Neumann conditions are imposed on the 
temperature field, which was not the case in our simula- 
tions. However, a parameter based on Q instead of AT 
is preferable in astrophysical applications because heat 
fluxes are better constrained by observations than ver- 
tical temperature differences. We will therefore seek a 
relation between Nu* and Raf*. Furthermore, in a flow 
dominated by rotation, which is necessarily nearly two 
dimensional, it seems plausible that heat flow through a 
plane z — const, should be determined solely by the dy- 
namics in that plane. Q would then be independent of 
the layer height d. If our working hypothesis is correct 
that Q is independent of n and v, and assuming Nu* is 
given by a power law, one has to find a scaling of the 
form Nu* oc Ra ^ . If in addition Q is independent of d, 
one has to find (3 — 1/2. 

Fig. 2] shows Nu* as a function of Raf*. The figure 
contains only those points for which 0.5 < Re Pr Ek 1 / 2 < 
10. This transition interval is small and does not corrob- 
orate any power law Nu* cx Ra^ at fixed Ek and Pr. 
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However, the data for different Ek and Pr collectively 
define an envelope which we regard to be the genuine 
scaling obeyed by the Nusselt number in the transition 
regime. The best fit to the data in fig. Q] yields 

Nu* = 0.11 • Ra° f f 5 (7) 

The exponent j3 = 0.55 ± 0.01 is measurably different 
from 1/2. There is some scatter in the points in fig. [4] 
around the power law ([7]). This scatter can be reduced 
by retaining data from a smaller interval of Re Pr Ek 1 ' 2 , 
so that the data are less affected by scalings valid in the 
neighboring intervals. 

Ref. [13| investigates thermal convection in a rotat- 
ing spherical shell. In this geometry, convection occurs 
mostly outside a cylinder tangent to the inner core and 
coaxial with the rotation axis, whereas the flow velocities 
are much smaller inside the tangent cylinder. Gravita- 
tional acceleration varies radially in the simulations in 
ref. [l3| and there is a zonal flow along circles of con- 
stant latitude which has no analog in our simulations. 
Despite all these differences, the heat flux in the spheri- 
cal geometry obeys Nu* = 0.077 • Ra J /* according to ref. 
[HI and the best fit to a compilation of data in ref. fl4T ] 
yields Nu* = 0.08 • Ra° f ^ 5 . The exponent in appears 
to be very robust. 

It is also interesting to draw a parallel with dimensional 
arguments for non-rotating convection [l5l |. If the heat 
transfer is independent of the layer thickness because it 
is determined by boundary layer dynamics, Nu has to 
behave like Nu oc Ra 1 / 3 . This exponent is generally 
not observed experimentally because of the presence of a 
large scale circulation. The assumption that heat trans- 
port is independent of thermal diffusivity and kinematic 
viscosity leads without rotation to Nu oc (Ra Pr ) 1 / 2 . 



While this scaling has been found in simulations avoiding 
boundary layers [l|| it remains elusive in any bounded 
geometry. In rotating convection, fig. 0] shows that a 
power law independent of diffusivities is a useful fit to 
the data, but the heat flow still depends on the layer 
depth. 

In summary, three different regimes of convection could 
be identified as a function of Re Pr Ek 1 / 2 . For small and 
large values of Re Pr Ek 1 / 2 , one approaches asymptoti- 
cally the scalings valid for rotating convection near onset 
and non-rotating convection, respectively. The cross-over 
occurs in a transition interval around Re Pr Ek 1 / 2 = 2. 
This contradicts the naive expectation that the transi- 
tion should occur when the Rossby number Ro = Re Ek 
equals 1. Even though Re is not a control parameter, 
the transition criterion is useful when observations yield 
some information about the flow velocities in a celestial 
body. A case in point is the Earth's core, for which mag- 
netic secular variations provide us with estimates of typ- 
ical flow velocities around 5 x 10~ 4 to/s. Together with 
£1 = 7.29 x 10 _5 s" 1 and the generally accepted mate- 
rial properties inside the core of k = 3 x 10 _6 m 2 /s and 
i/=5x 10" 7 m 2 /s [13, one finds RePrEk 1 / 2 w 5, which 
places the Earth's core inside the transition interval. If 
on the other hand the Earth's core is driven by compo- 
sitional convection, a diffusivity of 7 x 10~ 9 m 2 /s should 
be used [13], leading to RePrEk 1 / 2 » 6 x 10 3 . 
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